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ABSTRACT 

We predict the biasing and clustering properties of galaxy clusters that are expected to be ob- 
served in the catalogues produced by two forthcoming X-ray and Sunyaev-Zel'dovich effect 
surveys. We study a set of flat cosmological models where the primordial density probabil- 
ity distribution shows deviations from Gaussianity in agreement with current observational 
bounds form the background radiation. We consider both local and equilateral shapes for the 
primordial bispectrum in non-Gaussian models. The two catalogues investigated are those 
produced by the eROSITA wide survey and from a survey based on South Pole Telescope 
observations. It turns out that both the bias and observed power spectrum of galaxy clusters 
are severely affected in non-Gaussian models with local shape of the primordial bispectrum, 
especially at large scales. On the other hand, models with equilateral shape of the primordial 
bispectrum show only a mild effect at all scales, that is difficult to be detected with cluster- 
ing observations. Between the two catalogues, the one performing better is the eROSITA one, 
since it contains only the largest masses, that are more sensitive to primordial non-Gaussianity. 



1 INTRODUCTION 



According to the concordance scenario of structure formation, the 
virialized cosmic structures that we observe today are the result of 
gravitationally-induced growth of small, primordial density fluctu- 
ations. The mean amplitude of these seed fluctuations can be de- 
rived by the level of Cosmic Microwave Background (CMB) tem- 
perature anisotropies, and makes up one of the cases for the exis- 
tence of collisionless dark matter. The question of how these pri- 
mordial fluctuations are generated is generically answered by as- 
suming a phase of accelerated (inflationary) expansion of the Uni- 
verse right after the Big Bang, that amplified quantum fluctuations 
of the matter field to produce the seed fluctuations. 

Many different models of inflation have been proposed over 
the years (Guth 1981; Linde 1983), differing mainly in the number 
and properties of the scalar field(s) driving the accelerated expan- 
sion. A generic prediction of many of these models is that the power 
spectrum of the primordial density fluctuations should be scale- 
free, with a spectral index very close to the Harrison-Zel'dovich 
value of unity. Other models predict instead a spectral index de- 
pendence on the scale (running index). Also, not all models of in- 
flation predict that primordial density fluctuations are Gaussianly 
distributed , forecasting instead a probability distribut ion that is not 
Gaussian jBartolo et all |2004| : lUoVerde et al.l l2008 i). CMB con- 
straints on the level of non-Gaussianity do not yet rule out this 
possibility, hence it is natural to search for alternative probes of 
it and to explore the observable consequences that deviations from 
Gaussianity can have on the formation of cosmic structures. Ad- 
ditionally, it is possible for the non-Gaussian amplitude to depend 



on scale, hence exploring cluster and galaxy scales, that are very 
different from CMB ones, can help constraining this dependence. 

Generically, it is expected that a probability distribution for 
density fluctuations that is positively skewed would produce a 
larger number of high density peaks, hence leading to the for- 
mation of more numerous massive structures. Similarly, a nega- 
tively skewed distribution would have the opposite effect. The con- 
sequences of this on the mass function of galaxy clusters, that 
are naturally the most affected structures, can be substanti al, and 
have been studied in a series of works ( Matarrese et aL 200(]|: 
Verde et al.'2000VMathis et al.'2004VKanH et al. 2007: Grossi et al.l 
2007; Maggiore & Riotto 2009; Grossi et al. 2009). Another, more 
recently studied effect of primordial non-Gaussianity on non-linear 
structures is on their spatial distribution. It is expected that in pri- 
mordial density fluctuations fields with probability density different 
from Gaussian the density peaks are differently clustered together, 
leading to different biases and correlation functions of galaxies and 
galaxy clusters. Only recently has a coherent framework for com- 
puting the correction to linear bias in gene ric models with prim or- 
dial non-Gaussianity been developed ( Mat arrese & Verd3l2008l) . 

In the present paper we employ these results on cluster bias- 
ing and earlier results on the cluster mass function suitably cali- 
brated against fully numerical «-body simulations in order to com- 
pute the clustering properties of galaxy clusters as would be mea- 
sured in the cluster catalogues produced by two forthcoming sur- 
veys with various models and levels of non-Gaussianity. The cluster 
catalogues that we ad opted are two amon gst those that have been 
addressed recently by iFedeli et"al] ( 1200 8l) . and they are produced 
by eROSITA and the South Pole Telescope (SPT henceforth). The 
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former is a space-based X-ray observatory, while the latter is a tele- 
scope aimed at millimetric an d sub-millimetric observations of the 
ISunvaev & ZeldovichI l ll972h (SZ) thermal distortion of the CMB 
spectrum. Both are expected to represent the state of the art instru- 
mentation in the respective fields for the forthcoming decade. 

The rest of this work is organized as follows. In Section[2]we 
detail the models of primordial non-Gaussianity adopted, including 
the way in which we computed the mass function and the linear bias 
of galaxy clusters in those models. In Section[3]we briefly describe 
the two cluster catalogues investigated in this paper, specifying the 
parameters of the related surveys and the scaling relations between 
cluster mass and X-ray/SZ observables that have been employed. In 
Section |4] we present our main results, showing the effective bias, 
observed power spectrum and spatial correlation functions for the 
non-Gaussian models and comparing them with the Gaussian case. 
Finally, in Section[5]we summarize our conclusions. The reference 
Gaussian model we considered is a standard ACDM cosmology 
with best fit parameters taken from the 5-years WMAP data re- 
lease, in conjunction w ith type la supernovae and Baryon Acous - 
tic Oscillation datasets jPunklev et alj2009l : lKomatsu et al.ll2009h . 
The present values of the density parameters for matter, dark en- 
ergy and baryons are f2„, o = 0.279, Q^.o = 0.721 and £l(„o = 0.046 
respectively. The Hubble constant reads Hq = hWQ km s"' Mpc"', 
with h = 0.701. The normalization of the cold dark matter power 
spectrum of primordial density fluctuations is fixed by erg = 0.817, 
and the slope thereof is n = 0.96. 



2 PRIMORDIAL NON-GAUSSIANITY 

A particularly simple and useful way to parametrize primordial 
non-Gaussianity consists in writing the Bardeen's gauge invariant 
potential <I> as the sum of a linear Gaussian term and a non-linear 
sec ond-order term that encapsulates the deviation from Gaussian- 
ity dS alopek & Bond"l990|; lOangui et al.lll994 IVerde et aljlioool : 
iKoiiiatsu & Spergel 2001), 



(1) 



In Eq. ([TJ the dimensionless parameter /nl, that weights the 
quadratic correction to the Gaussian random field Oq, is in general 
scale and configuration dependent. The symbol * denotes standard 
convolution, and in the particular case in which /nl is constant, it 
reduces to simple multiplication. Bardeen's potential O, on scales 
smaller than the Hubble radius, equals minus the usual Newtonian 
gravitational potential. 



A s recently noted by di fferent authors jAfshordi & ToUev 

20081 : IPillepichetai] l2008l : ICarbone et al] l2008l : iGrossi et al 



20091) . there is some ambiguity in the normalization of Eq. (TJ. Ac- 
cording to the Large Scale Structure (LSS) convention, that is the 
one used here, O is linearly extrapolated at z = 0. In the CMB con- 
vention instead cD is primordial, so that /nl = 8(+°°)f^i^'^ / si^) - 
1.3/™^, where g(z) is the linear growth suppression factor for cos- 
mological models different from the Einstein-de Sitter one. It is 
defined by 



1 8(z) 
l+zg(0)' 



(2) 



where D+(z) is the linear growth factor. This means that any con- 
straint on the value of /nl gathered from CMB data should be in- 
creased by ~ 30% in order to comply with the convention adopted 
in this work. 



Embracing the nomenclature o f lMatarrese & Verdd feOOSh . let 
us write down the Fourier transform of the present-time linear over- 
density filtered on some physical scale R. 



SR(k) 



2 T(k)¥ 



iy«(A:)0(k) = M«(fe)<I)(k), 



(3) 



where Wnik) is the Fourier transform of the top-hat smoothing 
function, T(k) is t he matter transfer fun ction and k = ||k||. For T(k) 
we adopt the fit of lBardeen et al.l (Il9861 with the correction due to 
baryon physics reported in SugivamJll995h . More sophisticated 
fits for the effect of baryons on the matter power spectrum exist 
(Eisenstein & Hu 1998), however we checked that these additional 
refinements are unimportant at the scales of interest here. 

From the above Eq. lO, it follows that the relation between 
the power spectrum of matter density fluctuations extrapolated at 
present, P(k), and the power spectrum of the Newtonian potential, 
P^(k), reads 



P{k)THk) 



2 T(k)¥ 



3 H^Qn 



P<t(k). 



(4) 



As a consequence, if the primordial matter power spectrum is scale- 
free, P(k) = Ak" as in our case, the potential power spectrum can 
be rewritten as 



P.s,ik) : 



k"- 



Bk"- 



(5) 



Several models of inflation predict that the bispectrum of pri- 
mordial perturbatio ns in the potential ass umes a particular shape 
that is called local iLo Verde et Zll2008h . and it is such that the 
magnitude of the bispectrum itself is maximum when one of the 
three momenta (ki,k2,k3) has a much smaller magnitude than 
the other two ("squeezed" configuration). In these models /nl is 
a dimensionless const ant, and the bispectrum can be written as 
dCreminelh et al.ll2007h 

Bo(ki, k2, kj) = 2/nl52 [l^l-^kT" + + ^r'^r'] • (6) 

Other kinds of inflationary scenarios predict a primordial bis- 
pectrum with equilateral shape, in the sense that it is maximized 
by configurations where the three arguments have approximately 
the same magnitude. In the latter case, the primordial bispectrum 
takes the cumbersome form 
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(7) 



Most importantly, in inflationary models that predict an equilat- 
eral primordial bispectrum, the parameter /nl is in general depen- 
dent on the sca l es. We adopt here the functional form suggested by 
iLo Verde et ai]( l2008l) . according to which 



/NL(kl, k2, kj) = /nl.O 



k, + k:, + k-t 



(8) 



The functional form of Eq. ^ is chosen in order to avoid vio- 
lating the WMAP constraints, in the sense that /nl.o represents 
the non-linear parameter evaluated at the scale ^cmb = 0.086/? 
Mpc~' roughly corre sponding to the largest multipole used by 
iKomatsu et aiT ( l2009l) to estimate non-Gaussianity in the WMAP 
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Figure 1. The correction to the mass function for different kinds of non-Gaussian initial conditions. The top four panels refer to a primordial bispectrum with 
local shape, where the black solid lines have /nl = -12 and the red dot-dashed curves refer to /nl = 145. In the four bottom panels we show results for 
bispectra with equilateral shape, where /nl = -200 for the bottom group of lines and /ml = 330 for the topmost group. Different lines refer to different scale 
dependence for the non-Gaussian amplitude /nl: k = (black solid and red dot-short dashed lines), *: = -0.1 (blue dotted and magenta dot-long dashed lines) 
and *: = -0.2 (green dashed and cyan long-short dashed lines). The four leftmost panels show the amplitude of the correction as a function of mass for two 
fixed redshifts, as labelled, while the four rightmost ones show the con'ection as a function of redshift for two fixed masses. 



data, / = 700. The constant free parameter k is assumed to 
be \k\ <s: 1 between CMB and cluster sca les. Consistently with 



ILo Verde etall l l2008h : ICrociani et alj ilOO^ we assume small and 
negative values for k, that enhance non-Gaussianity on scales 
smaller than CMB. We adopted k = 0,-0.1, -0.2. 



2.1 Mass function 

Gene ralizations to non-Gaussian models of the standard 
IPress & Schechter (1974) mas s function hav e been presented 
in iMatarrese et al.l ilOOOt) and ILo Verde etal] ( |2008|) . Both ap- 
proaches ass ume that deviat i ons fr om Gaussianity are small. 
In particular, IMatarrese et al.l ( I2OO0I) use the saddle point ap- 
proximation to compute the probability distribution of threshold 
cross ing, and then truncate the resulting expression to the skew- 
ness. |Lo Verde et al ] (l2008h instead approximate the probability 
density function for the smoothed dark-matter density field using 
the Edgeworth expansion and then perform the integral of the 
probability distribution for threshold crossing exactly on the first 
few terms of the expansion itself. The two approaches give quite 
similar results, and both have been shown to give reasonable 
agreement with ful l numerical simulations of structure formation 
jGrossi et al]|2009h . provided the linear overdensity threshold for 
collapse is corrected for ellipsoidal density perturbations according 
to Ac — > Ac ^/q, with q = 0.75 (see also^aggiore & Riotto 200§). 
In this work, we ad opted the formula of | Lo Verde et alj j2008l) . 
following which the IPress & Schechteil h974h mass function for 
cosmologies with non-Gaussian initial conditions can be written 
as, setting S^iz) = ^c/D+(z), 
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(9) 



In Eq. 1(9} (Tm is the rms of the density fluctuations field 
smoothed on the scale R corresponding to mass M, while SsiM) 
is the normalized skewness of the same field. It reads SiiM) = 
~/NL,oy"3(M)/o"M' where in the local case obviously /nl.o = /nl^ 
and the third-order moment jiiiM) can be computed as 



y"3(M) 
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(10) 



Under the assumption that the non-Gaussian correction to the 
mass function is independent of the approach that is taken to evalu- 
ate the mass function itself, the structure abundance in a cosmology 
with non-Gaussian initial conditions can be computed in compli- 
ance to a generic prescription according to 



n{M,z) = n^°\M,z) 



, nps(M,z) 
nS,\M,z)' 



(11) 



In the previous Eq. dlU. n'j,JiM,z) is the mas s function com- 
puted according to the iPress & Schechteil U9l4l formula, while 
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vr \M,z) is the one computed as specified by t he preferred pre- 
scripti on, that in our case was the one detailed in lSheth & Tormeiil 
("2002), both of them evaluated within the Gaussian model. The 
[Press & Schechter ( 1974) mass function in the non-Gaussian cos- 
mologies, «ps(M, z), can be computed analytically following Eq. 

In Figure[T]we report the corrections to be applied to the Gaus- 
sian mass function in order to obtain the non-Gaussian one. As 
mentioned above, in the equilateral case we adopt three different 
values for the exponent k, namely k = Q, -0.1 and -0.2. As for the 
values of /nl in the local model and /nl.o in the equilateral mod- 
els, we adopt the co nstraints given by the 5-years WMAP dataset 
jKomatsu et al. ^2009). that are the tightest presently available, with 
the exception of the work of Smith et al.. C20p9,) that however quot e 
limits only for the local shape. According to Komatsu et alj ^20091) . 
for local non-Gaussianity is allowed to vary between -9 and 
111 on CMB scale. Since we are adopting the LSS convention, we 
used/NL = -12and/NL = 145 as extremal values. Similarly, for the 
equilateral cases the CMB constraints give -151 ^ /NL(f ^ 253, 
hence we adopted the extremal values /nl.o = -200 and /nl.o = 
330. 

Different theoretical studies, based both on analytic and nu- 
merical investigations, have addressed the capability of different 
observables in constraining the non-Gaussian amplitude, suc h as 
the a b undance of massive virialized structur e s (^ atarrese e t al.l 
2000'; Verde et al. 2000'; Mathis et al. 2004; Kang et al. 20q1 
Grossi et al. 2007) , halo biasing (Dalal et al. 2008; McDonald 20o5 
and this work), galaxy bispectrum Isefusatti & Komatsul 120071; 
Jeong & KomatsJ |2009|). mass density distribution tarossi et al.l 
20081) and topology jMatsubaral l2003l; iHi kage et al. 2008), inte 



_grated Sachs-Wolfe effect jAfshordi & Tollevll2008l;ICarbone et al 
20083 , Lyo- flux from low- density intergalactic medium (viel et al 



reionization jCrociani 



20091) . 21 -c m fluctuations (ICoorayll2006l ; |prileDich et al.ll2007h and 
' etalj|2009i) . 



In general, all these methods 
provide weaker constraints than the CMB data, with an interesting 
exception being Isiosar et alj ( |2008|) . 

As can be seen in Figure [T] the correction to the mass func- 
tion increases both with mass and redshift. At z = it can be up 
to ~ 25% for the most extreme masses, irrespective of the chosen 
shape for the primordial bispectrum. At z = 1 the corrections raise 
up to a factor of ~ 2.5 for the local shape and the equilateral one 
with positive /nl.o- For negative /nl.o, the corrections become arbi- 
trarily large for the most extreme masses and high redshifts. This 
fact is not a concern, since at very high masses and redshifts, the 
abundance of objects is virtually zero. 



2.2 Bias 

Primordial density fluctuations with non-Gaussian probability dis- 
tribution cause a scale-dependent modification to the linear bias 
for give n mass and redshift. Herea fter, we adopt the approach de- 
tailed in iMatarrese & Verdd | |200^ . according to which we define 
the function TR{k) as 



o-l Jo 

r' s*(<r, VS,fc) 



dfi 



(12) 



where ctr is the rms of density fluctuations filtered on the scale R, 
a = + k^ + 2(kfi and B^, is the bispectrum of the non-Gaussian 



potential, where the three arguments have been replaced by scalars, 
since S4>(ki, k2, ks) = B^(ki,k2,ki) for both non-Gaussian shapes 
considered in this work. An important feature to be explored of 
cosmological models with non-Gaussian initial conditions is the 
configuration dependence of higher order correlation functions. 

Given all the above, the Eulerian bias in models with non- 
Gaussian initial conditions can be written as 



b(M,z,k) = l + bf\M,z) 



1 + 



Ah(M, z, k) 
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(13) 



where the Lagrangian bias in the Gaussian model was assumed to 
take the form 
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.[Z)+(z)0-MF'' + [V^AePj' 

(seel Mo & Whitelll996l ; lsheth & Tormerilll999l ; ISheth et al.ll200ll) . 

The parameters are here set to p = 0.3 and a = q = 0.75. The 
correction inside the square brackets in Eq. l |13t is 

Ab{M,z,k) A, TRik) 



bf\M,z) DAz)MR{ky 



(15) 



In the particular case of a local primordial bispectrum, the re- 
lation for TRik) can be simplified to 



TR{k) 



8;r2o-: 
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C^MR{0P<t.iO X 



if 



(16) 



This prescription for the correction to the li near bias has been 
confro nted with «-body numerical simulations in lDesiacgues et al] 
( I2OO9I) , where it was found a disagreement between the theory and 
the numerical experiments for some ranges of bias and scale. How- 
ever, simila rly to what happen for the non-Gaussian mass function, 
ICrossi et"ai. (2009) found that instead a reasonable agreement can 
be reached with the positi on Aq — > Agi?, obtain ing results also in 
qualitative agreement with iPillepich et al] ( I2OO8I) . We adopted this 
position when computing the bias in non-Gaussian models in this 
work. 

In Figure|2]we show the correction to the linear bias as a func- 
tion of mass and scale. The redshift in this figure is always fixed 
at z = 0, since the correction scales simply as S^^iz). The shape 
of the correction factor as a function of the scal e is in qualitative 
agreement with the work of iTaruva et alj ( l2008t) . while a precise 
quantitative comparison cannot be made due to the different set of 
parameters that have been used. In particular, we note the expected 
fact that the correction to the li near bias in the case of local shape 
grows as oc k^^ at small k (Mat arrese & Ve rde 2008). Conversely, 
the correction decreases with decreasing scale in the case of equi- 
lateral shape. In the latter case, the correction is also much smaller, 
reaching at most ~ 10% at small scales and extreme masses. 



3 CLUSTER CATALOGUES 

Evaluating the clustering properties of galaxy clusters we referred 
to two forthcoming surveys, one in the X-ray band and the other 
one in the millimeter regime, exploiting the thermal SZ distortion. 
The first one is the wide eROSITA survey, while the second is the 
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Figure 2. The scale dependent correction to the hnear bias for different kinds of non-Gaussian initial conditions. The top four panels refer to a primordial 
bispectrum with local shape, where the black solid lines have /nl = -12 and the red dot-dashed curves refer to /nl = 145. In the four bottom panels we show 
results for bispectra with equilateral shape, where /nl = -200 for the bottom group of lines and /nl = 330 for the topmost group. Different lines refer to 
different scale dependence for the non-Gaussian amplitude: k = Q (black solid and red dot-short dashed lines), k = -0. 1 (blue dotted and magenta dot-long 
dashed lines) and k = -0.2 (green dashed and cyan short-long dashed lines). The four leftmost panels show the amplitude of the correction as a function of 
mass for two fixed scales, as labelled, while the four lightmost ones show the correction as a function of scale for two fixed masses. The redshift is always set 
to z = 0, since it can be easily scaled. 



SPT survey. These are the two most promising in order to distin- 
guish models with a strong redshift evolution of dark-energy by 
using the cluster co rrelation functions, out of the five considered in 
iFedeli et aI.U 2008Y 

The eROSITA wide survey is planned to have a sky cover- 
age of ~ 2 X lO'' square degrees down to a limiting X-ray flux 
of ^"lim = 3.3 X lO"''* erg s"' cm"^ in the energy band [0.5, 2.0] 
keV ( see also the dark-energy task force white paper lHaiman et all 
l2005t) . In order to convert this limiting flux into a minimum mass 
at fixed redshift we employ ed the s et of scaling relations de- 
scribed inlFedeli et alj ilOO^ (see also lSartelmann & Whitel2003l : 
iFedeli & Bartelmanij|2007ir They consist of the virial relation be- 
tween mass and X -ray temperature with norm alization based on 
the simulations by 'Mathi esen & Evr^S ( 1200 ll ). together with the 
luminosity-temperatur e relation required by the observations of 
lAllen&Fabi^ ( ll998h . These imply a relation between mass and 
bolometric X-ray luminosity of the kind 



L(M,z)- 



3.087 X lO-^^ergs^'r^ 



M 



lO'^Mf 



■h(z) 



(17) 



where the mass is expressed in units of Mq. In this deriva- 
tion it is implicitly assumed that the luminosity-temperature 
relation do es not evolve with redsh ift , as justified by th e 
studies of Mushotzkv & Scharj il99H) : iReichart et all ( Il999l) : 
iHashimoto et alj ilOOA and that the steepening of this relation 
at galaxy group scales d oes not apply ( Osmond & Ponman 2004 ; 
iKhosroshahi et alj2007h . 

We converted the bolometric luminosity implied by Eq. l |17| l 
i nto a band luminosi t y by m odeling the intra-cluster plasma with 
a iRavmond & Smitlil jl977h model with a metal abundance Z = 



0.3Zo jFukazawa etatll 19981 ; ISchindleilll999l) . The plasma model 
has b een implemented with the xspec software package 1 Arnau3 
The luminosity is then converted into a flux by using the 
luminosity distance in the appropriate cosmology. 

For the SPT catalogue the predicted survey area is of ~ 
4 X 10^ square degrees a nd we adopted the specifics detailed by 
iMaiumdar & MohrH2003h . according to which the limiting SZ flux 



density at a frequency vo = 150 GHz is Sy 



5 miy. It is 



likely that eventually the survey area will be larger than predicted 
in .Majumdar & Mohr (.2003.) (M. Joy, private communication, see 
also Staniszewski et al.[ 2008h , however this is relevant only for 
evaluating the uncertainties on the observed correlation functions. 

In order to link the minimum SZ flux density with a mini- 
mum catalogue mass, we used the scaling relation between mass 
and Compton ^j-parameter integrated over the solid angle cov ered 
by the virial sphere of the cluster given bv lSehgal et al. ( l2007h . to- 
gether with the relation between the integrated Compton parameter 
and the nominal SZ flux density. The result is 



S„o(M,z): 



M 



2.592 X 10**mJy 
(Da(z)/1Mpc)2 \ lO'^M, 



E(zy 



2/3 



(18) 



where the mass is again expressed in units of Mq, DpXz) is the angu- 
lar diameter distance out to redshift z and Eiz) = h(z)/h. We remark 
that not all the features of the scaling relations described above are 
well established, especially concerning their redshift evolution. We 
however believe that they are the most suitable given our aims. 

The minimum mass included in the eROSITA and SPT cat- 
alogues as a function of redshift for a variety of dark-energy 
cosmologies (incl uding standard ACDM model) is shown in 
iFedeli et al.l ( |2008|) . Since the scaling relations adopted here de- 
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Figure 3. Left panel. The all-sky equivalent redshift distribution of galaxy clusters in the catalogues produced by the eROSITA and SPT surveys. The Gaussian 
model is represented by the thin black line. Right panel. The ratio between the redshift distributions of clusters inside eROSITA and SPT catalogues obtained 
in different models with non-Gaussian initial conditions to the Gaussian case (thin black line). Line types and colors are the same as in previous Figures. 



pend only on the geometry of the Universe, the minimum mass in 
each of the non-Gaussian models adopted in this work is the same. 
In Figure [3] we show the redshift distributions for the eROSITA 
and SPT catalogues in the various models with non-Gaussian initial 
conditions, as well as the ratio thereof with respect to the standard 
ACDM cosmology. The (all-sky equivalent) redshift distribution is 
defined as 

Jr>M2 
n{M,z)dM, (19) 
Ml 

where n(M, z) is the differential mass function and G(z) is the vol- 
ume contained in the unit redshift, that in a flat Universe can be 
expressed as the Jacobian determinant 

G(z) = r-(z)^(z). (20) 
dz 

In Eq. J20b the function r(z) is the comoving radial distance out to 
redshift z- 

The difference between the two catalogues is evident, in that 
the redshift distribution drops to zero already at z < 2 in the 
eROSITA catalogue, while it is still significant at z ~ 3 in the SPT 
one. As already discussed bv lFedeli et al. (20081), this is due to the 
different redshift dependence of the scaling relations adopted, in 
particular by the fact that the X-ray flux drops as the square of the 
luminosity distance, while the SZ flux density drops as the square 
of the angular-diameter distance. 

The difference between different initial condition models are 
mostly visible in the right panel of Figure[3] showing the ratio with 
respect to the standard ACDM cosmology. In general, for both local 
and equilateral shapes of the primordial bispectrum, differences are 
more enhanced in the eROSITA catalogue than in the SPT one, 
and this is an obvious consequence of the fact that the minimum 
mass included in the former is larger at any given redshift. Since the 
deviations from the Gaussian mass function increase with mass (see 
Figure[TJ, it is expected that the corresponding redshift distribution 
is more sensitive. 



For the case of local bispectrum with negative /nl> the devi- 
ations with respect to the Gaussian model are always very small, 
due to the fact that in this case /nl is very close to zero. On the 
other hand, the deviations are more appreciable when /nl is pos- 
itive. For the equilateral shape and /nl < 0, the departures from 
the Gaussian cosmology become arbitrarily large in the eROSITA 
catalogue, however this happens at z > 2, where the number of ob- 
jects in the catalogue is practically vanishing. If we limit analysis at 
z < 2 for the eROSITA catalogue we can see that deviations from 
the Gaussian redshift distributions are at most of a factor of ~ 2.5 
in this catalogue. For the SPT catalogue, departures from the Gaus- 
sian redshift distribution reach up to ~ 80% at the highest redshifts, 
z ~ 3, where we still have < 10^ objects in the catalogue. 



4 RESULTS 

In this section we first of all assess the evolution of the effective 
bias, both as a function of redshift and of scale, in case the initial 
conditions are not Gaussian. The effective bias is basically given 
by the linear "monochromatic" bias weighted for the object abun- 
dance, and can be written as 

4;rG(z) C"- 

feetf(z,fc) = — 77-- b{M,z,k)n{M,z)dM. (21) 

N{Z) J Ml 

In Eq. Ml and M2 are the extrema of the mass interval that is 
encompassed by the catalogue at hand. In the realistic situations we 
are dealing with, M[ is the minimum mass of a certain catalogue at 
the given redshift jFedeli et al ]|2008l) . while formally M2 = -1-00. 

The fundamental difference between this formula and the 
equivalent one for Gaussian initial conditions is that here the effec- 
tive bias has a dependence on the scale in addition to the redshift 
dependence. Such scale dependence is inherited by the correction 
to the monochromatic bias. It is interesting to explore how big this 
scale dependence is, in order to understand if it could be detectable 
by looking at the spatial distribution of clusters. In Figure |4] we 
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Figure 4. The effective bias of galaxy clusters included in the catalogues produced with the eROSITA and SPT surveys considered in the present work. The 
thin black lines refer to the Gaussian model, while the other lines refer to initial conditions of non-Gaussian type, with colors and line types equal to those 
used in previous Figures. The scale is fixed to k = 0.2h Mpc"' in the top group of panels and to k = ().05/i Mpc"' in the bottom group, as labeled. The four 
rightmost panels in both rows show the ratio between the non-Gaussian effective biases and the Gaussian one. 



show the redshift dependence of the effective bias for fixed scales, 
and the ratio thereof with respect to the Gaussian case. We selected 
the two fixed scales at k = 0.05/? Mpc"' and k = 0.2h Mpc"', in 
order to probe linear and mildly non-linear regimes. 

One first thing to note, that might seem counter-intuitive, is 
that those non-Gaussian models for which the abundance of objects 
is larger, i.e. those with positive /nl, are also those where galaxy 
clusters are more biased with respect to the underlying matter den- 
sity field. While this is an obvious consequence of the fact that the 
sign of the correction to the monochromatic bias given in Eq. l |15t 
depends on the sign of /nl, one would naively expect that in mod- 
els where it is easier for a density peak to collapse into a virialized 
structure, the structures would be less biased. This is in fact what 
we found in our previous dark-energy related work, iFedeli et"al] 



l l2008l) . However, with a little bit of attention, it turns out that this 
kind of reasoning is not correct. When the probability distribution 
for density fluctuations changes, the distribution of density peak 
heights is also changed, meaning that there are more (less) high 
peaks if /nl is positive (negative). Hence, the different abundance 
of structures reflects the different distribution of peak heights, not 
the threshold for structure collapse. At the same time, for positive 
/nl the density peaks themselves are also more clustered together 
with respect to the Gaussian case, implying a larger effective bias 
as found. 

Similarly to what happens for the redshift distributions, the 
differences between Gaussian and non-Gaussian models are more 
evident in the eROSITA catalogue as compared to the SPT cata- 
logue. This is a consequence of the fact that more massive objects 
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Figure 5. The effective bias of galaxy clusters included in the catalogues produced with the eROSITA and SPT surveys adopted in this work, relative to the 
scale-independent Gaussian effective bias. The thin black lines refer to the Gaussian model, while the other lines refer to initial conditions of non-Gaussian 
type, with colors and line style as in previous Figures. The redshift is fixed at z = (top-left group of panels), z = 0.5 (top-right), z = 1 (bottom left) and z = 2 
(bottom-right). 



are included in the former, that are more affected by non-Gaussian 
initial conditions. The increment due to non-Gaussianity can be up 
to ~ 20% for the effective bias. When the scale at which the effec- 
tive bias is evaluated is increased, the difference with respect to the 
Gaussian scenario generally decreases for the equilateral shape, in 
agreement with the general behavior to the bias correction exam- 
ined before. Likewise, we find a slight increase of the difference in 
the effective bias for the local models. 

In Figure[5]we show the ratio of the effective bias with respect 
to the Gaussian model for the two cluster catalogues at hand as a 
function of scale for different fixed redshifts, ranging from z = to 
z = 2. 

Reflecting the trend that has already been noticed in the cor- 
rection to the linear "monochromatic" bias, the ratio between the 



effective bias in a non-Gaussian model with local-shape bispectrum 
and the Gaussian one grows indefinitely at small values of k for pos- 
itive values of /nl, while it decreases below unity for negative /nl 
at the same scales. For instance, at z = and for /nl = 145, the non- 
Gaussian bias is already a factor of 2 larger than the Gaussian one 
at scales > 400 Mpc. At z = 2 the scale where this happens re- 
duces to ~ 100/7"' Mpc. The effect is milder for negative /ml since 
in this case the value of the non-Gaussian amplitude is substan- 
tially closer to zero. Analogously, the ratio between non-Gaussian 
effective bias computed in models with bispectra of the equilateral 
shape and the Gaussian case have an opposite trend compared to 
the local shape, in that it increases with the wavenumber k. Also, 
it has a much milder variation with respect to the local case, be- 
ing almost constant over the range of scales considered here. In the 
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Figure 6. The observed power spectrum for galaxy clusters inside the eROSITA and SPT catalogues, as labelled in the plots. The two leftmost panels show 
results for non-Gaussian models with local shape, while the rightmost two show results for models with equilateral shape relative to the Gaussian case. In all 
panels, the thin black lines refer to the Gaussian model, while other line types and colors follow the same convention as in the previous Figures. 



equilateral case, the deviation with respect to the Gaussian case is 
only of ~ 1% at z = 0, and grows up to ~ 15% at z = 2. We also 
note the usual difference between catalogues, with the eROSITA 
one displaying larger variations than the SPT one. 

One can use the effective bias of galaxy clusters and the power 
spectrum of the dark matter density field to construct the power 
spectrum of clusters that should be observed in t he dif ferent sur- 
veys. Following the notation of 'Matarre se et al.l ( 1 19971) (see also 
iMoscardini et"5]|l998l I2OO0I [2OOI, 200^ we may write the ap- 
proximate expression 

PoUk) = -j -^bl^{k,z)P{k,z)dz, (22) 
where the normalization constant T reads 

In the two previous equations zi and zi are the limiting red- 
shifts of the cluster catalogue at hand. Practically, we shall have 
Zi - 0, while Z2 is the maximum redshift at which objects 
are pr esent in the catalogue. We adopted the IPeacock & Doddsl 
( Il994h fit for computing the nonlinear matter power spectrum, 
as we believe it suffices to our purposes (see the discussion in 
iFedeli et alj2003) . Additionally, we decided to neglect the r edshift- 
space distortion to be applied to the matter power spectrum jKaiseil 
1987';'Zaroubi & Hoffman" 199^;'Matsubara'2000'). As explained in 
Moscardini et al. ( 2000), th is co rrection results in a small change to 
the correlation function jBorgani et al.iri999i) , corresponding to an 
at most ~ 6% increase on the observed correlation length for deep 
surveys. We safely ignore this correction because it is not very de- 
pendent on the non-Gaussian model, and we are mainly interested 
in relative differences. 

In Figure |6] we show the results on the observed power spec- 
trum for the galaxy clusters contained in both the eROSITA and 
SPT catalogues, as well as for all the non-Gaussian cosmologies 
considered in this work. In that Figure, for non-Gaussianity of the 



equilateral shape we only show the ratio of the observed power 
spectrum to the Gaussian case, in order to better highlight the dif- 
ferences that would hardly be visible otherwise. Let us first focus on 
the local shape model, that is perhaps the most interesting one. On 
small scales, the non-Gaussian power spectrum always approaches 
the Gaussian one. Perfect coincidence is never achieved, since the 
corre ction to the linear bias never vanishes ^Matarrese & VercQ 
I2OO8I) , and even if it would, the redshift distributions would still be 
different. On large scales, the power spectrum for the non-Gaussian 
models with positive /nl increases without bound, according to the 
behavior of the correction to the "monochromatic" bias discussed 
above. Similarly, when /nl is negative the observed power spec- 
trum decreases far below the Gaussian one. 

In the non-Gaussian model with equilateral shape of the pri- 
mordial bispectrum, the differences with respect to the Gaussian 
model are very small, such that they are almost not visible unless 
we take a zoom of some region or we perform the ratio to the Gaus- 
sian case itself. This is in agreement with the behavior of the red- 
shift distribution and effective bias discussed above, and shows that 
these kinds of models should be more difficult to be distinguished 
by the Gaussian scenario using the spatial distribution of galaxy 
clusters. We additionally note that, coherently with the discussion 
presented in the previous sections, the effect of non-Gaussianity of 
equilateral shape on the observed cluster power spectrum is more 
marked for the eROSITA catalogue than for SPT. 

By performing the Fourier transform of Eq. M2\ with respect 
to wavenumber we then computed the observed correlation func- 
tion of galaxy clusters, that we denote with ^obs('')- In Figure|7]we 
report the ratio of this function to the Gaussian case for the various 
non-Gaussian models considered in this work and the two clus- 
ter catalogues we adopted. Errors on the correlation functions are 
computed via the bootstrap method, and are then propagated to the 
ratios. The first thing to note in this Figure is that the relative errors 
for the eROSITA catalogue are extremely small, and much smaller 
than those for the SPT catalogue. This is in part due to the dif- 
ferent area of the sky that is covered by the two surveys, with the 
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Figure 7. The ratio between the observed congelation functions for the various non-Gaussian models to the Gaussian case, represented by the thin black line, 
as a function of spatial separation. The top two panels refer to non-Gaussian models with local shape, while the other panels to models with equilateral shape, 
with different scale dependence of /nl, as labelled. The shaded regions denote errors computed according to the bootstrap method for each correlation function 
and then propagated to the ratios. Colors and line styles are the same as in previous Figures. 



eROSITA one being ~ 5 times larger than the SPT one. However, as 
noticed above, the SPT survey area might be underestimated here, 
and it is possible that eventually the errorbars for this catalogue will 
be smaller than depicted in Figure|7] 

For non-Gaussian models with primordial bispectrum of the 
local shape, the deviation of the observed correlation function with 
respect to the Gaussian case grows with increasing radius, a trend 
that reflects the one already observed in the power spectrum. With 
SPT, only the model with positive /nl can be reliably distinguished 
by the Gaussian case, with the errorbars being too large to allow 
the same for the model with negative /nl. On the other hand, with 
eROSITA the errorbars are small enough to allow the separation 
also of the latter models, if sufficiently large scales are probed. 
Considering the equilateral shape instead, the differences between 
Gaussian and non-Gaussian models are much more reduced, and 
maximized at scales < lOA"' Mpc. As a consequence, only in the 
eROSITA sample the models could be distinguished, while in the 
SPT one the errorbars would be too large. An exception to this is 



maybe given by the model with positive /nl,o and the most extreme 
scale dependence of /nl, namely k = -0.2. In this case the devia- 
tion between models might just be large enough to be resolved. 

A popular way to rapidly quantify the correlation strength is 
the correlation length tq, defined such that ^obs('"o) = 1- The ef- 
fect of non-Gaussianity with positive amplitude /nl results in an 
increase in the measured correlation length of ~ 20% - 30% with 
respect to the Gaussian case, with the precise value depending on 
the model and on the catalogue considered. 

It is interesting to note that the ratio of non-Gaussian correla- 
tion functions to the Gaussian one is almost constant for the equi- 
lateral shape and for the local shape at spatial separations < 20/j"' 
Mpc. At these scales, the difference with respect to the Gaussian 
case might be accounted for by a change in the normalization erg 
of the primordial power spectrum, that enters quadratically in the 
normalization of the observed correlation function. Turning the ar- 
gument around, this gives the precision with which is necessary to 
know (Tg in order to disentangle the effect of non-Gaussianity. It 
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Figure 8. The ratio of the spatial congelation functions observed in the different catalogues to the Gaussian case. Only clusters included in the redshift bins 
labeled in the plots are considered in each panel. For eROSITA we consider only two redshift bins, while for SPT we consider three of them, as detailed in 
the text. The thin black line refers to the Gaussian initial conditions, while other line types and color codings are the same as in previous Figures. Errors are 
computed with the bootstrap method for each coirelation function and then propagated to the ratios. 



turns out that for models with positive /nl this precision is of the 
order of ~ 10% - 15%, depending on the model, a precision that is 
already available. 

As a final step, we computed the observed power spectrum 
and spatial correlation functions when only clusters within selected 
redshift bins are considered for each catalogue. For the eROSITA 
survey we considered two bins, z < 0. 1 and z > 0.1, while for 
the SPT survey we adopt the three bins, z < 0.1, 0.1 < z < 0.3 
and z > 0.3, thanks to its far wider redshift distribution of objects. 
As explained in Fedeli et al. ( 2008), this choice insures an approx- 
imately equal number of pairs of objects in each bin. In Figure [8] 
we show the ratio of the correlation functions obtained for the dif- 
ferent kinds of non-Gaussianity assumed in the present work to the 
Gaussian case. Let us focus first on the eROSITA catalogue. In the 



low-z bin we basically have no signal, since the difference between 
non-Gaussian and Gaussian models is very small and the errors are 
quite large. Instead, for the redshift bin z > 0. 1 the deviations from 
the Gaussian case are large and the relative errors are small, allow- 
ing a significant separation. It should be noted that in this latter 
redshift bin, the absolute errors are actually slightly larger than in 
the former, however the correlation function is also larger, so that 
the relative error is effectively smaller. These conclusions apply to 
both shapes of the primordial non-Gaussian bispectrum. 

Focusing on the SPT catalogue, we note that the deviations of 
non-Gaussian models from the Gaussian one increase with increas- 
ing redshift, and the size of the relative errors decreases accord- 
ingly. The non-Gaussian models with equilateral shape can never 
be distinguished from the Gaussian case, since the deviations there- 
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from are always too small compared to the errors. The situation is 
different for the local shape, where deviations from Gaussianity can 
be detected in the highest redshift bin (z > 0.3) and at sufficiently 
large spatial scales. The latter however only applies to the case with 
positive /nl, since the other one is still too similar to the Gaussian 
model. 



5 SUMMARY AND CONCLUSIONS 

In this work we evaluated the main properties of galaxy cluster cat- 
alogues to be built with the two forthcoming survey performed with 
eROSITA and SPT, exploring cosmological models with various 
kinds of non-Gaussian initial conditions. In particular, we focused 
on the observable redshift distributions, on the effective bias and on 
the observed power spectrum of galaxy clusters obtained with the 
two catalogues in the different cosmologies. The two cluster cata- 
logues are constructed adopting the predicted survey properties and 
simple yet realistic scaling relations between mass and X-ray/SZ 
observables. 

The non-Gaussian models adopted have both local and equi- 
lateral shape of the primordial bispectrum, with different amplitude 
of the non-Gaussian deviation, compatible with the bounds coming 
from CMB and other probes. The redshift distribution of objects 
in the two catalogues is only mildly affected by primordial non- 
Gaussianity, resulting only in at most ~ 80% modification for SPT 
and a factor of ~ 2.5 for eROSITA at the most extreme masses and 
redshifts. The reason for this difference in the two catalogues is 
that X-ray flux drops more steeply than SZ flux density with red- 
shift, hence the latter catalogue is more dominated by high-mass 
objects compared to the former, which in turn are more affected by 
non-Gaussianity. 

The effective bias is affected in a way coherent with previ- 
ous work. Namely, it displays a scale dependence that is absent 
in Gaussian models. The deviation of the effective bias with re- 
spect to the Gaussian case grows at large scales in models with 
local shape, while it mildly grows at intermediate-small scales for 
models with equilateral shape. As a consequence, while the devia- 
tions from the Gaussian case can be very large in models with lo- 
cal non-Gaussianity if the wavenumber is small enough, in case of 
non-Gaussianity with equilateral shape maximal deviations range 
from ~ 1% at z = up to ~ 15% at high redshift. It is interest- 
ing to note that non-Gaussian models with a positive /nl provide a 
larger abundance of massive structures and a larger effective bias, 
meaning that not only in these models the large peaks that even- 
tually collapse into bound structures are more numerous, but that 
also peaks themselves are more clustered together. The opposite 
obviously applies to the case of negative /nl- 

The power spectrum that is predicted to be observed with the 
use of the two cluster catalogues above reflects the behavior of the 
effective bias. As a matter of fact, for non-Gaussian models with 
equilateral shape, the power spectrum is very similar to the Gaus- 
sian one, with maximal deviations occurring at intermediate scales 
and reaching up to ~ 20% for the model with k = -0.2. On the 
contrary, the cluster power spectrum deviates significantly from the 
Gaussian case for non-Gaussianity with local shape. In particular, 
when /nl is positive, the power spectrum grows indefinitely at large 
scales, being already ~ 2 orders of magnitude larger than the Gaus- 
sian power spectrum at A: ~ 3 x 10"^^ /i Mpc"' . For the same kind of 
models but negative /ml, the power spectrum decreases far below 
the Gaussian one. All of these conclusions apply quite indepen- 
dently of the catalogue adopted, except that in the eROSITA one 



the effect of non-Gaussianity tends to be slightly more marked than 
for the SPT catalogue, as a consequence of the different mass com- 
position of the two, as explained above. 

By computing the Fourier transform of the observed power 
spectrum we also evaluated the observed correlation function that 
is expected to be measured with the catalogues detailed above in 
the different non-Gaussian models. We estimated the expected er- 
rors on the observed correlation function by using the bootstrap 
method. Coherently with the behavior of the power spectra we find 
that the observed correlation functions for the local non-Gaussian 
models deviates strongly from the Gaussian model at large spatial 
separations. The deviation is expectedly more marked for the model 
with positive amplitude of the non-Gaussian contribution, since in 
that case /nl is more distant from zero. For non-Gaussianity of 
equilateral shape instead, the deviations with respect to the Gaus- 
sian case stay always quite limited, and never grow above ~ 20%. 
The relative errors on the observed correlation function for the 
eROSITA catalogue are much smaller than those for the SPT cat- 
alogue, mainly due to the largest area of the sky that the for- 
mer cover. Therefore, not only the differences between models are 
slightly more enhanced in eROSITA compared to SPT, but also the 
errors are smaller in the former. This certainly makes eROSITA the 
ideal tool for this kind of study. 

We also demonstrated that subdividing the two cluster cata- 
logues in different redshift bins, the better results are always got 
when the highest bins are considered. This is consistent with devi- 
ations from Gaussianity being larger at higher redshift and higher 
masses (that are preferentially selected at high z). Even in this way 
however, only eROSITA seems to be able to detect deviations from 
Gaussian initial conditions. Separating a given cluster catalogue in 
different redshift bins does not give any particular advantage over 
considering the entire sample, however it demonstrates the impor- 
tance of including high-z objects, that produce the bulk of the sig- 
nal. 

Two possible distinctive signatures of non-Gaussian initial 
conditions that have not been discussed here are somewhat re- 
lated to the present work. The first is the statistics of voids in the 
large scale structure . This issue has been recently add ressed by 
iKamionkowski et"^ ( |2009|) (see also lGrossi et 'Z||2008^, showing 
that the abundance of large empty regions can indeed be used to put 
constraints on the non-Gaussian amplitude /nl at the level of few 
tens. The second is the use of maximum cluster mass as a func- 
tion of redshift as a discriminator between models. While this is 
an interesting issue that deserves exploration, the comparison with 
real observations would be quite difficult, since objects at the ex- 
treme mass end are especially rare. Also, because of this paucity, 
the scaling relations at such high masses are not well defined. 

Before concluding, it is worth mentioning that the errors esti- 
mated in this work do not take into account the presence of scatter 
around the scaling relations used to link the dark-matter halo mass 
with cluster observables. This scatter will have the effect of some- 
what increase the size of errorbars, so that our reported values are 
likely to be lower limits. Still, in order to produce large deviations 
from our results, the scatter would need to be extremely skewed, 
and more complete datasets would be needed in order to under- 
stand whether this is indeed the case. 

A complete statistical analysis of the predictive power of 
forthcoming cluster surveys in terms of shape of the primordial bis- 
pectrum and level of non-Gaussianity is certainly a step to perform, 
but goes beyond the main purpose of this work. We can conclude 
that the effect of primordial non-Gaussianity on the clustering prop- 
erties of galaxy clusters is genetically mild, but depends strongly 
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on the shape of the primordial bispectrum that is chosen. It is likely 
that constraints competitive with those from the CMB can be given 
in this way using the eROSITA catalogue. For the SPT catalogue 
this is probably not possible, unless the survey area is increased 
such that the error on the observed spectrum can be reduced below 
~ 10%, or clustering measurements can be pushed out to very large 
scales, but in this case, only if the non-Gaussianity is effectively of 
the local shape. 
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